Self-Consistent Nuclear Shell-Model Calculation Starting from a Realistic NN 
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First self-consistent realistic shell-model calculation for the light p-shell nuclei is performed, start- 
ing from the high-precision nucleon-nucleon (NN) CD-Bonn potential. This realistic potential is 
renormalized deriving a low-momentum AW potential V| ow _k that preserves exactly the two-nucleon 
low-energy physics. This Vi ow -k is suitable to derive a self-consistent Hartree-Fock basis that is em- 
ployed to derive both effective single-particle energies and residual two-body matrix elements for the 
shell-model hamiltonian. Results obtained show the reliability of such a fundamental microscopic 
approach. 

PACS numbers: 21.60.Gs, 21.30.Fe, 27.20.+n 



The nuclear shell model is the foundation upon which 
the understanding of the main features of the structure 
of the atomic nucleus is based. In such a frame, a cen- 
tral role is performed by the auxiliary one-body potential 
U, which has to be introduced in order to break up the 
nuclear hamiltonian as the sum of a one-body compo- 
nent Ho, which describes the independent motion of the 
nucleons, and a residual interaction Hi: 



v 2 

H = J2 V H = T+v = (T+U)+{V-U) = H +Hi 



z<3 



(1) 

Once Ho has been introduced, it is possible to define a re- 
duced Hilbert space, the model space, in terms of a finite 
subset of ffo's eigenvectors. In this way the unfeasible 
task of diagonalizing the many-body hamiltonian (1) in 
a infinite Hilbert space, may be reduced to the one of 
solving an eigenvalue problem for an effective hamilto- 
nian in a finite model space. 

During the last forty years, a lot of efforts have been de- 
voted to derive effective shell-model hamiltonians, start- 
ing from free nucleon-nucleon (NN) potentials Vnn, 
which reproduce with extreme accuracy the NN scat- 
tering data and the deuteron properties (see for instance 
0). In particular, during the last ten years many re- 
alistic shell-model calculations have been performed to 
describe quantitatively with a great success theproper- 
ties of nuclei over a wide mass range 0, 03, IH E| 

In all these studies, for the sake of simplicity, a har- 
monic oscillator potential has been adopted as auxil- 
iary potential. This choice simplifies the computation 
of the two-body interaction matrix elements, as well as 
the derivation of the effective hamiltonian, by way of a 
degenerate time-dependent pertubation theory 0, |g] • 

A more fundamental microscopic choice for U is the 
Hartree-Fock (HF) potential, that is self-consistently 
derived from the free potential Vnn- Such a choice 
leads to a so-called self-consistent shell model |9|. It 



is well known, however, that modern realistic NN po- 
tentials cannot be used directly to calculate a HF self- 
consistent potential, owing to their strong repulsion at 
short distances. A traditional approach to this problem 
is the Brueckner-Hartree-Fock (BHF) procedure, where 
the self-consistent potential is defined in terms of the re- 
action matrix vertices G |10|. However, this approach 
cannot be considered the basis for a fully self-consistent 
realistic shell-model calculation [jj, because the choice 
of the BHF potential for states above the Fermi surface 
cannot be uniquely defined [TT| . 

■ Recently, a new technique to renormalize the short- 
range behavior of a realistic NN potential by integrating 
out its high-momentum components has been introduced 
|12|. The resulting low-momentum potential Vi ow _k is a 
smooth potential that preserves the low-energy physics of 
Vnn and can be used directly to derive a self-consistent 
HF potential This paves the way to perform a full 
self-consistent realistic shell-model calculation. 

In this Letter, we present results of such a kind of 
shell-model calculation for light p-shell nuclei, starting 
from the high precision CD-Bonn NN potential [T^ |. 




FIG. 1: Behavior of -Ehf with hui and N for He. 

First, for the sake of clarity, we outline the path 
through which our calculation winds up. The first step 
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consists in deriving from the CD-Bonn potential a low- 
momentum potential Vj ow _k defined within a cutoff mo- 
mentum A. This is done in the spirit of the renormaliza- 
tion group theory so obtaining a smooth potential which 
preserves exactly the on-shell properties of the original 
Vnn E2- In Ref. a detailed discussion has been done 
about the value of A and a criterion for its choice. Ac- 
cording to this criterion, we have used here A — 2.1 fm _1 . 
The so obtained Vj ow -k, hermitized using the procedure 
based on Cholcski decomposition suggested in Ref. [l5| . 
has been used directly to solve the HF equations for 4 He 
doubly closed-shell core. We remove the spurious center- 
of-mass kinetic energy [Tt| writing the kinetic energy op- 
erator T as 



tonian is written in operator form as 



r = icL5>-Pi) a ■ 



(2) 



So, the hamiltonian can be re-written as 



H 



PI 



A 

E 2^ 



Pi ' Pj 

mA 



(3) 



The details of our HF procedure may be found in 
Ref. H3- I n our calculations the HF single-particle (SP) 
states are expanded in a finite series of N = 3 harmonic- 
oscillator wave-functions for 4 He. This truncation is suf- 
ficient to ensure that the HF results do not significantly 
depend on the variation of the oscillator constant fav, as 
shown in Fig. 1, where the behavior of the HF ground- 
state energy versus fku has been reported for different 
values of N. The value of fao adopted here is 18 MeV, 
as derived from the expression fao = 45A~ 1 / 3 — 25A~ 2 / 3 

m 



TABLE I: Experimental [2^1 and theoretical SP energies (re- 
ferred to 4 He closed-shell core) calculated including up to 
second- and third-order diagrams in VW-k- The experimen- 
tal widths of the states are also given. Energies are in MeV. 





2nd order 


3rd order 


Expt. 


r 




1.109 


1.355 


1.97 


1.23 


nOp 1/2 


4.800 


4.970 


3.46 


6.60 




-0.081 


0.202 


0.89 


0.65 


vOPl/2 


3.796 


3.984 


2.16 


5.57 



The HF SP eigenvectors define the basis to be used for 
our shell-model calculation. More precisely, we assume 
that doubly-magic 4 He is an inert core and let the valence 
nucleons occupy the two HF orbits 0p3/2 and Opl/2. 

The next step is to derive the effective hamiltonian 
for the chosen model space. Starting from the time- 
dependent perturbation theory 0, the effective hamil- 



H cS = Q-Q> / Q+Qi / Q / Q-Q> 



Q Q+ .. 



(4) 

where Q is the irreducible vertex function Q-box com- 
posed of irreducible valence-linked diagrams, and the in- 
tegral sign represents a generalized folding operation. Q' 
is obtained from Q by removing terms of first order in 

Vlow-k- 



TABLE II: Experimental and calculated energy levels, in 
MeV, in 6 Li and 6 He. The calculated results obtained in- 
cluding diagrams up to second and third order in V[ ow _k are 
compared. See text for comments. 



6 Li 


Expt. 


3rd order 


2nd order 


Binding energy 


31.995 


31.501 


32.198 


£ x (l+0) 


0.0 


0.0 


0.0 


£ x (3 + 0) 


2.186 


1.912 


1.757 


£ x (0+l) 


3.563 


2.627 


2.401 


£ x (2+0) 


4.312 


5.215 


5.170 


£x(2+l) 


5.366 


4.911 


4.590 


£ x (l+0) 


5.65 


7.825 


7.305 


6 He 


Expt. 


3rd order 


2nd order 


Binding energy 


29.269 


29.797 


30.837 


£x(0+l) 


0.0 


0.0 


0.0 


£x(2+l) 


1.8 


2.108 


2.072 



We take the Q-box to be composed of one- and two- 
body diagrams through third order in Vi ow -k- 

After the Q-box is calculated, the energy-independent 
H e g is then obtained by summing up the folded-diagram 
series of Eq. (4) to all orders using the Lee-Suzuki itera- 
tion method y|. However, such a method is appropriate 
for degenerate model spaces only, this is the case when 
using a harmonic oscillator basis, as pointed out before. 

The model space we are dealing with is non-degenerate, 
the unperturbed HF SP energies (respect to 4 He closed- 
shell core) being 4.217 and 7.500 MeV for 7r0p 3 / 2 and 
nOpi/2, 3.106 and 6.577 MeV for v0p 3 / 2 and vOpi/2, re- 
spectively. So we use a generalization of the Lee-Suzuki 
iteration method to sum up the folded diagram series, 
expressed, in the case of non-degenerate model spaces, in 
terms of the multi-energy Q-boxes |l8j . To our knowl- 
edge this is the first time that the above technique has 
been employed to derive a realistic shell-model effective 
interaction. 

As mentioned before, H c g contains one-body contri- 
butions, the sum of all these contributions (the so-called 
5-box 19]]) is what actually determines the SP energies 
that should be used in a shell-model calculation, ft is 
customary, however, to use a subtraction procedure [l9| 
so that only the two body terms of H D g are retained and 
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TABLE III: Experimental and calculated energies (MeV) and 
electromagnetic properties in 7 Li. The reduced transition 
probabilities are expressed in W.u.. 



7 I i 




v_^aic 


Binding energy 




oy.ooD 




n n 


n n 
u.u 


f / 1 - M 




U.4yo 


^(21 2) 




o.yyo 


M21 2> 


D.DU4 


0.99o 


-M22 2) 


/ .404 


7 01 q 


-^"422 2> 


0. 1 u 


o.yoo 


^22 2 > 


Q DQ 


Q 8QQ 


-Ml 2 D 


9.57 


9.564 


£ *(fi 1) 


11.24 


9.206 


Qgs [e mb] 


-40.6(8) 


-24.4 


A*gs H 


+3.256 


+4.28 




19.7(1.2) 


9.71 


B(E2;2- ^ §-) 


4.2 


4.18 


B(M1;|- - §;) 


2.75(14) 


2.37 



the SP energies are taken from the experimental data 
regarding the low lying spectra of the nuclei with one 
neutron or proton outside the inert core. In this calcu- 
lation we have followed a more fundamental approach, 
where the SP energies are the theoretical ones obtained 
from the S-box calculation (see Table I). 

TABLE IV: Experimental and calculated electromagnetic 
properties in 6 Li. The reduced transition probabilities are 
expressed in W.u.. 



b Li Expt. Calc. 

Q gs [e mb] -0.818(17) -0.442 

^i gs [nm] +0.822 +0.866 

B(E2;3+0 -»• 1+0) 16.5(1.3) 6.61 

B(E2;2+0 1+0) 6.8(3.5) 6.45 

B(M1;0+1 -> 1+0) 8.62(18) 9.01 

B(M1;2+1 -> 1+0) 0.083(15) 0.154 



In Tables II- VIII we compare experimental binding en- 
ergies, low-energy spectra, and electromagnetic proper- 
ties of 6 Li, 6 He, 7 Li, 7 Be, 8 Be, 8 Li, and 8 B USEHH 
with calculated ones. All calculations have been per- 
formed using the OXBASH shell-model code j^. A 
quantitative amount of data for different nuclei is taken 
into account in order to verify the reliability of a self- 
consistent shell-model calculation. As regards binding 
energy results, it is well known that shell-model calcula- 
tions give ground state energies referred to the closed- 
shell core. In our approach, we can consistently cal- 
culate the 4 He binding energy by means of the Gold- 
stone linked-cluster expansion |13| . So, our theoretical 



TABLE V: Same as Table IV, but for 7 Be. 



Be 


Expt. 


Calc. 


Binding energy 


37.600 


37.751 


-M21 2) 


0.0 


0.0 


E-x ( h ^ h } 
x V 2 1 2 / 


0.429 


0.469 


-Mil I) 


4.57(5) 


3.922 


§) 


6.73(10) 


5.782 


& (f 2 I) 


7.21(6) 


7.123 


Mgs [nm] 


-1.398(15) 


-1.013 


B(M1;£- - 1") 


2.07(27) 


1.81 



binding energies are obtained summing the ground state 
energies of the open-shell nuclei to the calculated 4 He 
binding energy, whose value is 25.967 MeV including up 
to third-order contributions. Electromagnetic properties 
have been calculated using effective operators |2J] which 
take into account core-polarization effects. 



TABLE VI: Experimental and calculated energies (MeV) in 
8 Be. 



8 Be 


Expt. 


Calc. 


Binding energy 
S x (0+0) 
E42+0) 
£ x (4+0) 
S x (2+1) 
i? x (2+0) 
£x(l + l) 
£ x (l+0) 
£ x (3+l) 
S x (3+0) 
i? x (4+0) 


56.50 
0.0 
3.04 
11.40 
16.63 
16.92 
17.64 
18.15 
19.01 
19.24 
19.86 


54.010 

0.0 
2.985 
10.056 
15.889 
15.174 
15.927 
15.070 
18.863 
16.983 
21.191 


TABLE VII: 


Same as Table IV, but for 8 Li. 




8 Li 


Expt. 


Calc. 


Binding energy 
£x(2+l) 

£x(l + l) 

£ x (3+l) 
E x (l+1) 
Sx(4+1) 
£x(0+2) 


41.276 

0.0 
0.981 
2.255 
3.21 
6.53 

10.822 


44.365 

0.0 
0.922 
2.074 
4.667 
6.929 
8.230 


Q gs [e mb] 
/i g s [nm] 
B(M1;1+1 -> 2+1) 
B(M1;3+1 -> 2+1) 


24(2) 
+1.653 
2.8(9) 
0.29(13) 


26.7 
+2.89 
2.67 
0.38 



In Table II we compare calculated binding energies and 
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spectra of 6 Li and 6 He, obtained including contributions 
up to second- and third-order in perturbation theory. In 
all other tables calculated quantities refer to a third-order 
H cS . 

From the inspection of Tables II- VIII, we can con- 
clude that the overall agreement between theory and 
experiment may be considered to be quite satisfactory. 
This agreement is of the same quality, and in some cases 
even better, of that obtained in our previous work [25| . 
where a realistic shell-model calculation for p-shell nu- 
clei was carried out within the framework of the semi- 
phenomenological two-frequency shell model. It is worth 
to note that the inclusion of a realistic three-body force 
could lead to an overall improvement in the agreement 
with experiment, as shown in Refs. [2H l27j 



TABLE VIII: Experimental and calculated energies (MeV) 
and electromagnetic properties in 8 B. 



8 B 


Expt. 


Calc. 


Binding energy 


37.74 


40.429 


£ x (2+l) 


0.0 


0.0 


£ x (3+l) 


2.32 


2.061 


£ x (0+2) 


10.619 


8.245 


Q gs [e mb] 


64.6(1.5) 


40.1 


Mgs H 


+1.036 


+1.73 
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